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ABSTRACT 

Tbe uwchAiusia tbiaugh vlucb SMtef-a 2 ««i bouldets giov to b&-ei 2 ed ]»ifl"a>jcni»>h 
m protoplAiutary discs is s eobjeci of Active resMreb, SDce it e aiUcai for plsnat 
foraksbon. lb Awid spir&lljag tbe protostu due to AerodysAink drag, objects 
must rapidly grow from aa>sfted pebhlea, wtucb are tightly cooplad to tbe gas, to 
lar^ booldas of l-lOOm is diameter. It is already wall knowo that over-dea^ks io 
the gaasoQS of tbs «<m<* provide potential stes for the coUectioo of solids, 

oiwi that sigu&raot density structures in tiu gaseoew of ths (S 4 ., epiraJ 

dessi^ wewa) can ti^ s^ids e£c«iitly enou^ foe the siriid CMopooeat of the disc 
to ondergo further grantatiooaJ collapse doe to tb^ own setf^gravity to work, 
we emp^y the eanciL ccob to conduct local flwing sheet simulaCoos of massive sef^ 
gravitating proloplanetary dkos, to stody the effect of aotacydosic Dana^it lortices, 
or eddies, on the evohitioo of eobds in theee discs. We find that these types of structures 
ere extieoaly «•»«"» ccecemrating grnftH iaiennediate-sised dust partwUs 
with fricbon timee comparable to, or leas than, the focal ortntal period cd the dkc. This 
lead to significaat {T^wr•den 8 ilies m the sobd compoosot of the disc, with dosty 
anhaiw a mant a comparable lo, 4 od €760 highsT, *!»»" thoee within epiraJ density wave^ 
meteasmg the rate of gravitatioiaal collapse of solids into boond sQnctures. 

Key wordar accrebon, accretion • gravitabon • bydrodynaoucs - lOEtabibUes - 
planets and satalUtss; formation 


1 INTRODUCTION 

leige-«cale spiral density waves cOBuonly arise In a peelcv 
planetary disc due to the gravltauoeal Inetablluha at the 
early eiage cd disc Bhs, vbaa U H in&tenLly maasve It baa 
been known tor some that the density wavts can act 
to UMii^cet >!■» oolWarda, alfowtng tasss 

to acoete onto the proloelar, and could weD t« the prlnusy 
transport provided thus by s elhgariey dut- 

S ieeachstetaMctfetKfomattoe (eg,|^.h^a^ 

IRiee etaPBma . 

The other pertubnUon t^e in that las rectived 
Bweb attendan, partly In eonneettoe wAb pi^--* friltiial ttc> 
an voUas They can arise in both self>grnvltaLlng 
and tuawelhgmvRaLlng discs by several naans. A < 

radial Itructun a the disc can trigger the R^by wave 
Inste^ty, to tarn nailts in the fafiB«iiA«< of ver- 

f ee.. ILovelace eteh^lU^aeaiLnnetal 
. Im55h et aL 12011 llid BMM IftsBlv etal 


tins 


Vortlee can aho be generaled as a result of other 
hydrod ymtnic s jst s 

Miu» Jkui.. b ■HM Staaen et al 1 1555^ 

n P»»uir.i—.tenni IR—oln et allapiat or tbe Kelvin. 


assss 


• p^gOroe ec ak 


Bgleihola kastaMUly ^ seed vorLlcIty stil^ 

B5^ Although the dynamics and evelnttoe cf wtlcsu In 
non>selPgravKaling dises has been unds active lnvtBtJga> 
tlon kd a wbils, ^^e gng idea o f their televanee to 
formatnn bvlB arae ilWflt . sysiemaUc 

atudlns ^ vwto dyaaciiles b s elfrg avtt allng cUsa are quite 
recent CfankssAvdl ii 103 IgSa Eivn 5^ l5ia ~!E3 
[201% « alllMlifl It wa 
cMBmcnly thought that spiral dsnslly waves are the only 
perttirbatUsi type fsTJrnt b selhgraltating dleca, and the 
mstence of other nwte aho participating in the overall 
disc dynamics had been l^cved In mwiy sti^he ol sdt* 
gnvHating discs 

Developments In the tbeoey of non-ssl^^avllatbg 
discs, howa^, showed >*»«■ wUm can be CUnmcly) eM» 
pled with, and esdte , spiral itguaHy w a ves due to the dlsch 
differential rotation ifaoda et all I300fr Johnstm fe Qammld 


S03RAS 




































2 Gi^botu et 


S5o5 P'iftalBMoJ 5558l TlUB «ftMl v«s 

^lUBMWaMd to bd not* m the t u—1> 

of salf>grwttv If gLa^^yfliPCOTi Mw» 

ovof » hsM 4Aal)*ait by |UMn£lAAvU!_|<,^BellibVty 
iMCffl 0 bo*«d tbM lA bet «<alleal node Is abo subjM lo 
tbe lii£maee U asU^avtby aad eafi eiblbll gnvtt^lo^ Ur* 
sta^ty tnsh bclora fooparable lo tbat d 

^BSity ■«*« 1\is ckac^ ki^lM tbat vertbas at? as In* 
ptftaU as density waves to a Mmplata dynamical 

{detura d aalPgnvKailng dba. 

Subteqnaal dynamla and awliilkm of 

vMtJes vafUs pralli in bob-caK>gravUMb^ and ealf* 
fia^latlBg eysteu In boci>«elf>gnvn^ln9 discs, onaD* 
«wj» (antUyebnle) vcHbeeB, soon after fomattai. sibm^^ 
amrg a MM aacb otb a and incwae m aaa lag. LiaLaL 
BoH IGodtm ti livKl Uga laiB SOabr b“Bod«>^MLaJ 

555 < ' 55* a- 

Bbene^3~ 5c5 Wben tba vertices ranb a few 
be^bt, ri'inifiifaislillllii bacooew ifoportant, eauamg 
to anlt sptral densKy waves • due to Lbe abevs-maot^Bed 
coii^lmg phaBomasoa • and, conssqu^Uy, alowly 

3 m ^ iijis_sjlii ^ r eeesal bi‘»'^*ed crtistsi nes» 

Bbss et aJbOOd . Podo tn J 
la spite d ibis sboefc dualpatfea. ntfeyc^ife voe> 
tfeas cao be consldefad loo^livwd. coboenl aunvtUM In 
dm Bowaves, iWy esn be subject to radial aJ^raLkm 

—''t/a tbe InstabU* 




tsv ilLaBif fe Panak>BoifD00ftl^^--< ■*< ■! i«13ft , tba lypl* 
c^ Un»e<sak d wbkb is lar^ «Im«h vortex tncncwar/a^al 
tine, as s^aO aa to daaunirtfen b y tbe nariftnmlatt''’^ In* 
ata^ty in magnetised regkma fi^a^^C ^H ISouh . 

Tbe vortan la restated uadat kaflueace 

of dfec asibgrssrtHf mn^i'— — wistfcaa ace favoured 


In aanoA ad*^avily acts to oppMs 1^ Inreraa ca» 
cade «1 spectral energy Ur laigar senfes. as octnrs In tbe 
ooa<sd&gtavlutlng cas, and, taataad, awttera il to anaDet 
acaks. Tie bcal sheet ainiiatfees cd radlelly in* 

rtrmtnred <1 n. In tb e »s—»ee d baro'~'t>»w or p^ tiy wave 
insta^lks) bv IblMaausabtill fe Elal fioBft denon* 
sitate that voRlms In a riirasi TTrsi]' nerttoturfeuleBce an 
tsaaslanv sbort<bvedftrncturss ludargolng repeating cycles 
cd festnalton, ^wtb to sttas coaqianble to local Jeans 
feagtik and rvrMiial md dastruetkB due to tbe efe 

feels ef sd*^avlly and Kaplarlaa sbaac TUs psecui leais a 
few orbital parnds, aod ranks In a vary dlflerant, tstaparad 
to ti»e a 0 CNself>gravitating ewe, avohilkwary plcuire, 'akb 
-«»»T anall Im ceganked verlkes (addles) in rartousatagas 
cd erolntlob. ral^ than tbe relstlv^ larger scale *aD ce> 
gsnlsed vwtiees giadnally grerwlng v ia slow mmgna O n tba 
otbet bai^ in global gtHinUsa^a 1:y ILvra et alTifeJOtfl ; £Jij 
I2013b6 an really stencinnd ai^ ibay 

[Mlo U-Lirm (wkbeni gnvttoiurbnfeoce), devciep vortken 
via Ressby wave knatablllly at vortaaaUy Hiinr^tun Dm lo 
wlkgnvlly. tbcsic voAks have a wtoMitbal extent 

()e., a*iwinii«i anvenu^sa) oapand 'aitb sloB* 

lar sHtlhMl a^gravtly What la neleworiby. bow e v ar , 
vDAlcs la both local aad global g<rnniia>«ia a^M 4 eon* 
mon fsBiure - they {serdtica ovar-densitlei at tbelr cenLres 
aiib tbe of Tbonia’a paianwter Q. Hen 

over-denkles an Inpoasd m tbe density vastninns due lo 
of deaaily 'anves enKted by tbtse vortices 


Xo naaove discs, lbe gmvitatkaal lastaldbilM have 
been abown to wuse the diac lo fragnmnl. po ssibly ksd» 
tng to tb e <tnTt ffSfiaTtrn of ifw giaat pi»*i»«e fem UflBft 
Rke el all GoM hodi Afe 

though a ujBBlve {wotoplabeiary dIac Is tbcugbt to ba 
pcesenl rt early tines in tbe alar fecnaiica procsn, II 
Is nol dm stub can eod fesi eaougb dna to 
radlaUve to directly farm pi^rte rka this 

fe»ftt'.>llw*lBolev etalJIMOi 
iftice fe Amiiaad P0C8) . Bowevsr. even If elective cool* 
bg Hr**** pneeni within dlsa are loo kng lo alfew gl* 
ant planets lo fem vis bagmotamn. as ma u tn u e d above. 

dfecs in vew 2‘**b ! aeltgrorttMlng b ihdr 
early stage* flLin fe PHanyltty, (iwS. scneHnes even rt 
iWT Tbun fPw.^r gt al]l2CltiaiAnd. e*« t> 

| 2 p 0 ^ In this ewe. Indirect rnrtnt of plai^ . 4 ^, 

of dust partkks In gwaous denrSty utrurturw 
■>ii as a result of gravltattonal knstabtlUy - could 
be at work As outlined above, such euiiettifes can be ellber 
due U> denaity waves ce to ariiaiu. ■w iK e s 

In low nm Sfbere selkgrartly is negllifble. uw« 
IkSB have bean a^an to have a gtguA-^i qp 

tM bml evolution of dust paiilefc a p^ >h«aanel^ 


ms 


GodS 5c5^TB alspba b~..| pfirtf [ td^euletalj 
iLvta fc Ltd Big . IPu et al 1 1201 * Bbu at aUlMldl . In tbwe 
slnnlatboB, not bdudlng se 4 f>|M«ky. duet ipnln with • 
aide range of abts s^ere aboern to Im (rapped wuhln tbe 
ceatew of tbe vertex suucturss with lbs Krengest ccercen* 
IraLlm oerurilng partbeku with felctloa/ftoppifig tinw 
ceoparabfe to tbe ceMUl pesbd In other words, a 
anmoth. auAdeatly bngdived vortax is liidoed abk w efc c * 
lively trap dust p—eww 14 m « 4 ra, pcesibly acceferadng 
1 iiMiatlOB. 

lb date, boawvet, tbe evduUon of dust parlkfes in adt* 
gravitating discs hw been pctmastly studied in a large scale 
setting, eoncatntlng on ^ elect spiral denlly waves, b 
rtiirnus conpesren i of tbe dlse, have on dtisl parlkks. 
lltke e* all laJOADOOfll used ^bal to show tb^ 

jT,iffiiT d spiral ^nslty waves In t be disc esn have a 
diamlK eftoet 00 tbe partkle layw. Ingka^^ 

Lbe dust partldw are narimri to in te ract wKb tbe diac 
sok^ rka an aetodynank drag fwce. Tbe paAldes, esp^ 
dally tbeve with stopping »!■'»'— cooparabk to tbe fecal ^ 
bUal parfed. are alwn to >»«*•■»>»* trapped wubb lbe denaily 
wave Mcuctursstbal tom in lbe diR due lo self* ^avkty TUa 
esnsea Lbe fecal denaliy of sobd aateela] wttbln Lbe densUy 
laav cs to Has by an aider of magnJlike n nne cenpand 
lothevssage gas pi > ni ir u gt^ul Lea 
Hive lo negative aenms tbe dessliy ivave sLrurtuie. erwung 
sufeKepkrkan vsfecitks cn one cf the wave, nd super* 
Ksplerian cn Lbe otbn As a resuk, tbe drag feme causa 
diM grains lo drift toward tbe deslLy/presaare ^ 

end of densILy waves. 

Q^_ei^ OOOefI generalised tbass by 

mg the ipssrkat^al effscla of tbe sobd parilcku tbemelves. 
wbkb alfew throe partlck enbascenents to beccene even 
mceefwewunced. rrsiiltlrglfi grrilathnal^binutd ffiimpt 
or plsoeteetaBala, (asnlng ta tbe dke More reesMiy. t brse 
rwiiRs 'We ^p— on ta lQlbboos et aLI ifeatfi u^ 
kng tdgbei resolution fa'^dlnai^^al (2D) fecal ibmmg 
amt attriiWvin. of gd^g tvlferf*^ SMiiuis dken, treated 

[StiRice et alllSinL 


kn related fecal atndks 


(^9003 HAS. UNRASOOO.(S*7? 


























































































{rspptn^ vorticei in ielf^^mwiaUng discs 3 


IPaMdaktoeS ISI 3 > . Mlb -i«r..a.>«l jiip«r>p*md«B 
te iwmiai dwfcpaftlde dynainla. T^m siudKs tsuAd that 

vidra] ^«ult}' w»*«i ar« uery afiosfil at Lnpplag 
tlvdy "t^«n b«iB tevml c«nLliBel«c up lo tii~iT~r iwH pw> 
tkitB, wboaa iloppkBg Ubhb att conpaiabto lo ^ loeal ct> 
bhalpuHcd Vary Ufldly batudclunealp*rtKl«BCM (arm 
lb (be creau cd (bse ipiral «bmb due M (be Mr<Bg ^av* 
UMlOBel etUbcLlao Leiacui pwf^ia at stKh cobcM* 
uttuaa Hen prucMM eao mtilt b a eubMamial popula* 
UOB o( lar$e, pfninitwIiTiel rtrni soUd ob|eets bbag t^ildly 
forswd lb (be eaily stagat cd dUe evelutden. pevlag (be aay 
fiM (etnainal asd gae gtabt pi«nt via tbe (r^^ 

>a*—i tc a e a ecjwtei nedel <« IMIm^ e< ■LlUmft , 

D(M partKlm cu abo beco ai e (rapped by prenura 
Mgrirlfil~d ariib (bBUcyclobk) vortleea abd, w nur- 
Uoaed above, tble prseen bae bees oiebeively lb 

oo&<ed^invua(lag He r«k tt varttcea b (rapplbg 
duel lb self'graviiMbg le kn ubdarMocd abd aeoda 
fiir(ber InveellgaUeB, becaun cd tbe dlfieraal (LaiacUr cd 
mtexevolnueb Mlb seli’giavILy. PaKlck tiy 

eMtkee, ee well m nrlea dynaala Iteelf b (be pm et^ 
<d aeU^avtlv, wn fleet addicnml by via 

global diet Hoaievar, b eebup v^tka an 

geoBated bf (be ftcnby wave Imablllly lb as olbt*wlee 
diac bad b«.*e a iees>lmd regulae aatun - aa lb (be 
oQB<ad&gnvlta(lBg eaae* l^ugb aeliiiulliBlebedue 

to Ibe eflecl cd ael&gnvUy. Here, we a alw 

iietloo - a eelf'ipV'Uatlbg disc b a quMl-etaady giavtto((ir> 
buleeletale, wKbou( CMtbMua dttvlbgby gloM 
v Hasaby wase bMabfllUea. aad aaalyee bebavkui cd 
^bedded par(kkb ueibg Ibe abeaflbg abeet appRKk 
We wsuU Ute to bcte ci'»i a ebar(s (ocnparabk 
to tbe diet bogbl M kee) Kngtti irnVe juai eueb a 
greetedurbuleai stale, ratbar tbaa a > ^i»~/cQbefm ebe, 
M infm^ et—w. ei a eK»y>w* •* fli ea (e g, iGaBMiBOll 

E cream el all IMWTl Ae alMwb by maftiatee 5 *i ^^^I« 
iMOdt. Ibw'iute CMtUli cd tnaaleat sbvw lived anall-aule 
vMtkes peoduebg over-^bee regloba ovabld ao debbly 
■rave iLruclum (a edbet words, «« ertaad UiB study by Uk 
cludkbg duel partkiai wlib a goal of uuderMndlBg ibe role 
tbeee vortical euuciurea play b dust evohilUA. D(M dynam* 
le« lb tbe praitaee cri oaly dnelty wavee In tbe abeartag aii—» 
bn beea o aplored lb our prevloua papers IQIbbobe at all 
gOlg POld l. where larger deenab ibm were ueei Here, we 
>-fa» a —»niiw gat« rkat nabke us (o toom lb m 

enk vertbal sLructum abdexasdae duel ira^ 
fiHig edibility agaluet (be backdrop of deaelty vauattow 
due lo (be deaaliy waves feboefcs) Tbk will li^ us to ew 
tablieb tbe rok ^ voeUen lb plabct ttemuoe process lb 
aslAgravitatlbg dHes. 

Tbs paper le ergaaksed ae bCowa. la SeeUoa 2, we oufc 
Uae Ibe physical sewup and Initial adopted in 

our dbulatioee. la Secilon fl, we (he eedultob of 

tbe duel parlleks b l^ (he vw(«t sunetuiee b (he gas. 
Finally, eununary abd dlecwsiobe are $vefi b 4 . 


bg diec ealf>gradty, *'>g»>!»»>’ wHh (ha eqn^lone of taouon 
of parllcka coupled (o tbe gee ihrougb eerodynande 
diag fwce We ako Include Ibe self>pavUy of paitlelm uec' 
eawry to tbelr coUepoe properties Ae a bu> 

UHSl^ (ooL we enpley Lbe ?SNca. com Q ll Is a shttb 
ord er fdal I blri o eda temporal flblte dlSerance rode 
feee IBtand enhurd fMOJI fee full details 1 Tbe rtwen. coea 
imt e aoM e ss aumencal ewpeispartkiss ^^abeM,«l^ 

B5B1^ 

la Lbe shearing sheet apfwodmatkai, cUet dyoainka Is 
studied b (be local CarieslM cocedlnate from cebtied al 
wTtmM gdimlal radhie, rg, ftoo the nntral object and ro* 
latbg with (be disc's frequeacy, fl, at Lble radius. 

In ihb frame with (be z-and y-axle tapoelively, 

b tbe r^el aad aelmuLbal dlrotlkms, tbe disc's dlflerenUal 
rotetlOD aanlfaste llself as ao aUarutbalperalU flow ebaiaC' 
leflaed by s IbearaMar, 9, of bMkgroond vebdLy akeg tbe 
s—asde, ua = (0, —^hs). T^ egtaUbriuiB surlve 
of gn, He, abd are af—«t«Hy unlfwm Since 

Ua dke IS cool, abd tberafore »*«»» the w^ort ratb k smell. 
U/ro < 1, whm H s e,/(i Is (be disc leak belgb, nd c, 
w tbe gas sound speed. iMal ibearbg sheet model is 
based on tbe si^AMUe of Use bi^ 2D taydrodybamK equw 
Uoos of raolkb to (be lowest order lb this small pajamet^ 
—u nlng that Use dwe Is ako rasor ihb (see e 4 ., [Gammiej 

iaomb 

Our alHoilatlni gpans (hS l e g t ttl —Lt/i < X & 

£,/2, -Ly/2 € p < by /2 As a eustcenary, we ed^ 
l^ Maadard sbeer-pebodle boundary well am 

peebated lo ilmiibg a^et al (e g.. [Gammy 

Bdmt h v~» 5 3 [gg^ . Naosdy, foe any vertebk 

tHfiniiiig ‘rkrPy with background flow nik 

IrMted, we have 

/{*. p, t) = /(* + tr, p - «Oi,l, 0, (* - beubdery) 

/(t. p, l) = /C*, p +<), (p - boundary) 

Tbe shear parameter Is g s L6 for tbe KfU«>i.«i rslaltob 
pro&k edoi^od le (hU paper 


8.1 Gas doBsIty 

In lbs local (be CMllailty equation for Use vsrtkally 
btegretsJ gas density S Is 

^ + V (Eu>-fft*^-/oCE)=0 (1) 

where u = <i«, ng) Is tbe gw volodly calelWe to the brnk 
ground Kepbslan shear flew ug Dtm to (be hlgb>mdsr bu* 
HwiHiai fli-atw the rewciL oodb U ako a dkf* 

flislob (eta, /d, lo eosun bumerical sleblUly ami cjqilurc 


/o=<oC^'E+rb<o VE) 

Here (he qunLliy (0 Is the ibock ittifciMi eocfidsbl < 


a MODEL 

Tb tarwesugme (be dynamics (d parufles in 

a solAgravn^lng protepkbetary dke, we solve (be 2D Im 
cal sbeertng s^et oquatttm (sr gas on a find ^id. kndud* 


CD = 0.^(a^I{-r-u)4l){ii:)* 


Sm Wlp / foods 
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4 Gibbons et of. 


vW« is a MotflaoL, dufwUrlalAg t he stwa<»fc e t atoc t 
^HfciflLiti B outllaed la Ap^aadU B <i ILvTa eTail iOOOft , 
aad St B LtIN, U nl) tte 


a.a Gas vabelty 

Tbe ««[ua(bb of tu (be gas reladve to lbs uoper> 

Uirbed Kt^kctaa Sow fW tbe fano 


^ 4 (u 7)u - s ~^r ~ ^ * B4 ^Oor^ 

E, u-v. 

■T* — 


.-V*+tiu), (2) 


■We P la Ibe mlkall}' iMe^ated ^oua, S* k tbe eur> 
deaslty ot particks, ii la (be graailatttaal poteattal pr» 
dnced Jobill}' by tbe gas and pertlde nirface deuttW faee 
6) llie ledi band aide «i equatWb ( 2 ) deaslbea (be 
advectlcea by tbe >ek>d(y Seld, u. Itself aad by tbe ueBA 
KepleHab Go*. Oo (be rtsbt side, t^ br8( term le 
tbe ^wua ftece, (be aseoad aad thtrd teru i mu hit i ii (be 
teree aad (be etfee( cd abaai, reapectlvdy, aad (be 
fiwiMb t«RB daserlW tbe badi>reK(iea surfed cd tbe gas 
by ( be dasl paf**‘**wi dite to ueody*'«*^k drag iMee (see 
e g., ILwa et all l2D0ft li abauea et l20Hft , Tbis ttfce is 
Cnperuoul to betereafi (be uWC' 

Uy «1 w aad (be gu velodty u aad ■ lavecaely 

{wefpcrtbeal to (be stopping, es GlctkD tdiu, ff.et parU> 
dan Tbe Gflb term upmuuti (be Sarae diM to self-gravny 
cd tbe syalsB Fbally, (be code n erpUeli vlsceslty 

tenn. C. 

f,=»^ 9 »u 4 jrr u 4 2S-nLE) 

♦ + S 4 rin • >4. 


3.4 Dust partiolas 


T^ dust partkltB a>e treated as e large ti au* 

BMSlcnl rafpeisparUeltB Udbansan et aLl2Cmfl Mtb irniWhiiii 
r, s CD (be gdd aed eelodtUn y, = (v^.V*) 

Wee to (^ unptMtubed Kepkrlan totatbu eelodty, Vp> = 
(0, -fAZa), «1 pankko IB tbe Iota] Garteetaa buoe Tbae 
are evolved aecMdlng to 




(d) 


dv. 

* 3 r 


= -20* X Vp 4 flf>v»,> - Vi> 


U-V. 


(B) 


Tbe Stsi (WD (ecina ^ (be rlgbl g|de of aqiistlcn (6) 
njiriurpt (be {We and (bsoOD-lnerllal torcedue to 

abear. Tbe (htrd tetto U tbe forte eurted oo (be parUdes 
due (0 tbe eunboa graWaUcAal poeenual ^ Tbe fourtb 
(erm dtKrlbeo (be drag {wee cDsrted by (^ gas oo (be pai* 
(Icka w fcif* ansae from (Im uoladiy dltferonce btsaass (be 
t»e. Uabbe (^ gm, (^ partiefee do not feel (be preenure 
fatee In (be co^, (be drag force oo tbe par llefesfrwB (begs 
M ralnilelrH ^ tbe gu vdodty Geld to tbe pc»> 

sUbo of (be partUk, using' 



Tbe back>rwttloa oo (be gat from (a egustlep (21 

eakulated by (be e dmiue outlbmd lA|J<bamQ,e(_^ iBOll 


3 b Sel^pavlty 

Tbe ^avUetional poteMkal lb (^ dynainkal equaUeoa (2) 
and <b) la eafetilaled by Isvorttaig Pedum eqnaUim far lA 
wbkb ">«uaina «a (be rlgbl ade (be 
surface denslUs lo a ibcs thin dise (eg.. 



■bleb coBtatau bo(b shear vlsecfa(y aod a bulk vfeemlty for 
resoMiig Eete 8 Is (be (racskas rato^frstt^ 




and G ^ tbocli vktoeMy analogouB to (be 

abc<b HiffnaU*. cneAdeBl ^ above, bo( wbb Dm 

replaced vm. 


2.3 Entropy 


Tbe ?SNC3L COM uaw eattopy, r, as Us (bensody* 
Bsme variable, latber tbaa Internal energy, V. T^ eqnnllM 
fot entropy evedntfoo Is 


Si' 


(u-D.-erta^ = 




2!b8'- 




+ /.W 


( 3 ) 

■We (be 8n( (am oo (be tlgU kde Is lbs viscous 
beating t erm sod (^ st n oW term ■ aa »*>f41cU 
falfawmalGsamflBoihilRieeeiaUlgilihif 


Slit. ■e iMiimt ibe eflecUve tp (o be eDn> 

•sot (hrougboul (be asd taM Us valne 

te s son*', «bkii (s suAofenily large far (be disc to avoid 
fratpoestaLloo and seule Into a quast-flendy self.regulMed 
•ete Tbe term on (be ngbt side, /gl»), is a a*!*-* 
disalptlW^ (erm aiaU^ia to oudlsed fas ite densUy. 


A* B 4«G(E - Ce 4 Er - Epo)4(s) 


(8) 


tbe Pw( Rnmer t^ **!^**^ ifci > metbe d nndUW In 
(^ suppleBeBtary materiel of U^iansen et aU QOgj]). Note 
(hst (be pestnrb^ gas, E - lo, and partkfe, E* — Epa, 
surface denaltlea equaUoo (G), aince qb^ (be ^sv* 
KaUooal fn-ugni.l assoclaied ailb (be perturbed modon 
(and bence daWty perturbetlCD) of boib (be gueous aad 
solid ecenpements d^ertnW gtevliy frece la equadou (2) 
end (S). Bere, (be Easily k Fourier tnnsfrMwd 

from W (e,y>’plBne to tbe (lu, bp)-plaoe oUbool the Inter* 
nmdla te ewordLaate (rarDfaramlbn perfmsmd by|Gamn^ 
iaOOlf . Pw (bu purpwe, a Meadard FPT dee bun 

adapted to allow far tbe Cut tba( Ibe r^lal weveBumlrer 
hr ct eneb ipaLlal Fourier dependa oo Uiu h 

Ap(|) a iic,(0) 4 fflbpl b orde r to »e*w^ t^ a—m-- g|i^ t 
bou^ary ccsrdltloas alanSiMcMSiaWli 1- Pi<dbOd81, 


3.8 Unite nW Inillnl Cooditlons 

We normalise (^ quanllLUs by setting r.a = B a £• a 1. 
Tbe (Ime and vekdty units are [t| a n~‘ and (o) a cm, m 
suhlng In tbe orbUal period, r«*a 2ff. Tbe imltodfenglb 
Is (be seek bslgbt, |It a fj a c,8/n. T^ Initial Toomre 
parametD Q a k equJ to L tbrougheut tbe 

'niB sets (ba ^evliadona] ceonient G a o~'. Tlw 
surface deWly ol gae k bUlel^ set to unity e v ei yw b ere b 
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Fi^oM 1. l^^BiUiiiiuc mtIms ^ aAa 4 D ahWt 

fuft cl Om pwada a SaodsM^ Tbe Ace 

1 m» cfccaly iitUud mW > gnicytt^)’ gtatiMtarbali tStte 
f—- kJ «ibelar^TMale deiiii^aa%« 4 ntMmaQsrKA(* 
uvuidsMC covcnuM An tB^^ulw t i ctica (cm 9 Tbe 
AftnboOos al dna ** mum Uae a chon i& Ftg Q 

tL« tli~n 'ns i4nw<«iK la A sqiiar« wUh dlmn* 

tfaac L, =L, = aad K dtvl^ IM a ^ c4 

fi,xN,s lias X 1024 e«lU. ‘noi c^ce of unlw sets (he 
Am-mu- leA^th Lt s 30 Ji/*Q •• 6 .S 7 ff Obd Ihe KeTtSfP«e»d> 
Ifigsrldedlslttisa^s (UM62A. 

The gtf vetodtr Said Is ifilUeDy pertubed b y eec e 
N»i«P tAbdon SucluiUeas vHh (he rm enplKude = 
10~^. We LAhe (he vlsecKly Hd caeAdents t« be 

k s 10~* Abd y,it s D,t = 10 0, IhB ensues Mrtr—4i^l 
ste^ty Kroas the slmh hoae, edlhoul wuhlbg (hem eul. 
Wi we LO* gpiu ewealy between Sve frKtloo tlB»u, 

f, s |0.0l,01,1,10, l6lJ|n~' The dust pertkles ere IbiLMy 
nbdomly pieced throu^bout (he srlik aero velodty 

Abd do Mt Intenct edlh (be £bS u other dnsi partkta for 
the flrsl 10 arhlsel perttds Al (hh pohit, obce (he gcs has 
aeulsd kmo a i^irntl TTffij stele, which is ibdepeBdaot of lie 
i»ini«i condltuu, (he dng fereas becwees ibediM pArtkha 
Abd gas are IbLrodueed The rtirm (a (hm evolved tor a 
Airtber 30 orhKA, be&Me the pailklea* seli’snvILy Is lbln> 
dnced. The value fix the duabto-gas surlecA debslty 
raUfi, C*e/Cs, Is (^ea lo he OJA fiar each paniefe specha 


a RESULTS 
3.1 Gae BvoliKioB 

The evohiUoB o( the gAeeoua ^ disc Is lb 

good r^rmrm with that oheetved lo prevttus aael^o*** 
Sndkw of aeh'gievUAdeg <*'•'“ b (be sheaibg sheet with 
m Inpoesd eeaetaat (kme The lolUel perturhAtkae 

grow Abd develop bio flucluatbu b vabdiy, sur> 

debstty Abd poMoUaL ^*-**« (hao firm, whkh proceed 
to ^at Ihe gas, while (he cooBbg acte to rednce (ha b(er> 
oal enetgy aad ebtroy^ ct the gas Denalt}' slruciuiws ^ 
vebp, which are sheared eiil by dlfiareaUal rolat^a. Once 
the system has completed seva^ aiblle, (be beaUDg geoar> 
Med hy the p bakbced by the cooUeg Urm, and (be 
sybem setdea late a fpil irnnfly eelf'tegtilAted state lb 


Flgura 2 . FV a«t/< r. n the state at the «mw 

tima M ta Fi^(T] MomarMB ■mafi.ecslc, trr^alaslr ■haped, siim 
^doalc %e)rtuM («d£es), mseepoadtag to otgmsve PV regions 
(UaA dow and arV dStae arms) Iwn develcped la the sheet 
The ofwMmse etnacnoes jwodaced I 9 these forUem are 
sta to the gaa detm^ held dtm is Fig Q] The C4l0u.m^ ia rw 
sUKt o d to the rwge — SD < f < 20 W hens «i^A^»/r(fUAhM 
(he pnsmM of the sorUol itncnrs 


LhB grAVUotuhnJsbt stale, the (hetmAl. klbettf abd gt e»- 
KaUobal ebergUa ^ Ihe disc are ob Avsage costaol over 
»!»"• T^ SAturatad value of the sheer suew paraBMder o 
w proporUmal to the b vete e -"«<Ug time, while the mb* 
plttude of Ihe -<^»»a< tv itactuathpae w (he squaie toat o< (h e 
bverse cooUag lOeedbe et ALII20t<t tElee el aJ EOII 
Figure IQehows (he sorlace desmuy of (he gm after 40 othlts 
b the fully developed quml>steAdy ^evltoturbuleBee. Bere, 
(be IsrgeMCale ahr^w dua to deasUy wavw are vkthle alcns 
with "<»iw scale eunetutes aasodaled wHh vQr(kas0 

(lybamks of th aaa vert^es b (be eavltotuib ukhi state 
was liiiiiellg'il'il by^attmashvfll,&,Ri^ BOOtfl . which we 
recap beJow 

'Ae vcrtkal stnmtnres b the gas are chAraetetleed by 
(be potetaial vceUdty (PV), 


f. KUj-> 


( 3 -g)Q 


-E- 

During (ha bust pbme, i^uiiJ arml* irnV ptdUlve and Begh> 
live P V ngtob* are strongly sbaared bte strips, but negative 
PV {aaUcycloBlc] regbas s(u( to wtep up Ibte varteX'Uke 
atwetuea dne l o (he Kelvln<HelfnholtA bstablluy 

QiUba^ iaP0?n . The positive PV tegbne remab sbearsJ 
bte stn^, showing no el^tf of vutax fwtoMkm dnang 
Lbe entire couas of evolutbn Only these antkyulMUc r» 
gkibs, hevlng negative PV, are able te survive b l^ how 
by teklng Lbe farm voKlcas, Lbough they are not as regw 
m (hme Ai«i«<»otn^ b >»*i"-eelf-gavl(atltg dlsa 
te£_JUui^s^Re£evj 1200* l l<«hKo«n f~ 


ISheaet Jl20Q6ll . These amalbr aca^. Irregular AntlrycfanU 
voetkes, w eddMa, daracterlmd by tegallvePV, which have 


^ Is the t^pmda, vw show, by calnlSnig the atBoeoRelarUn 
funrbm he the PV ddd ta the rpM prwlj note, that tte ihorw 
ml length far time eerdon w nmre than 10 tuaaa 

largir than the gBd mO eUe, biplying that they are enSdemly 
emoltwd b ^oraaulatMM 
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6 Gi^botu et 


In Ibe qnnMMt«»dv 9 t^ nre ckwl}' wen in (be 
PV SeU In FI;, g Tley a|^i«i vnHMkeB ( 9 (rt^ 

Uirw) vttb ft fl>tr<iM lo tbe g«s euflece ^eelt}' dne 
u nUpfirtflHIhj' aBd •elf>gmMv, M cfeeriy seen in Fl^ |T| 
One ennisy tliM la (be pcenci cdiel^^avlly, i^FVleld 
—•—< lu *inpcta«' Ifi tbe detally Tbb B In ciatnet 
te (be (nee Mtb muhlpte nntillaTarlj lotnrnnbg ■'Wlleo In 
tbe M^Mlf'^nvUeilBg ebenrlng abee(, wb«r« denelly 
■nwi (ebecto) abed b y tbem are seea lo tbe gM denaiy 
SeU l^aa » all BOOtfl . Canfarlai Fl^|T|aDdS «.« n^ 
(la( (baae ■‘^Ucse Mlb aBMUer, by abeolute anlue, PV 
£lve Ma l» blgber over>densuh9 >*»«•» wUb larges, 
atadule value. PV (Im eaaple, abHe r^lon la (be 
Mir 4 ar^ deoelty tnep near (be lower leA cesaer) Tbeee over- 
^***^»MM at — ^eat^e ate even larger (ban tbe denelly 
varlaU^ near (be shock frtata. 

la Lha qiijtf>ele»dv gravllo(iirb«Jenl Mate, (be vcr> 
tin nrrtifitaimt/iiiiilre<r slrucluren undagolag recnrring 
phases of twBMUoa, groatb to sBes cenperabk to a local 
Jeane and evesitUBJ sbearlng and dastruclkB due lo 
tbe akMs cf eelFgravlty (gravltal^nal loMaU* 

Uy) Mid ba^ground Keplerlao sb^t EmL phase ty|ikally 
lasts about tsio orbUal perUds or kes. As a reauU. lo self* 
^avUatlog dks, tbe cweraD dynaalral pletiiie of vartex 
evobuton la irregular coosisUng cJ mkag imtirlmi wrUcee 
« diSeitol evdnitonary stages asd, t^ebtre, wllb tarhaia 
sues up to i~^«i Jeans stele By cooirast, b tbe non* 
aalFgiavitallng easa, long*llved vestas ssrudures parsM fcs 
bundrsdscdorbKs via oHsglng ol aoalle vcstkes loto larger 
(&(■ uoUl eventually tbelr alae naebes tbe disc atak i<»ic*<i 
It «^«ti< be Botoi et*— tbe db: coolbg, vbleb we bave 
das crib cd ustog a sbnpk ctasiaat *"*Ji*»g {BMcdpUoa, 
b raabty Is due to r^atlve Icwes b om tbe dlse ^ 
psi^ an lis ^ysi csi propardes (eg, |Jcbnscg_lt,C^Bn^ 
Ji Obetoualy, effects disc ibsrcata 


rryiBHiSEE 


dyt^Ki. srbKb. b tura. tbrougb tbe barodblc sonree, ^ 
tenniiMs voetex (potenilal vwlklty) dynaBiles. Eoweots. In 
tbe presence of aelf'^’avlty tbe aUuallon Is scenee/bat dlf* 
fereot la tbe ease d a i^V'iiig uaw, despite bw 

Ing a approslBsatkA, tbe settlsn late (be quasi* 
aoady ^Mtoturboleu sate by a bcaJ bab 

aiwe between butleg aad eaeillng, wbkfa is qualSatlvely 
«s*uw to tXt wttb a»—vy realletic r^latlva (e.g. 

Lp.».K Q >ffy [SB PeJev M SI SSSi fs^s^rzrS 

[ffilljl % P ■ eiiidi tiiil (be evobuhs of PV b (be 
{Benence cdreelisk cooBng wUI else be a»-»<i»> (o tbn( wUb 
a i^liKg (One u Ittig as tbe disc rvsrdM b (be 

quutsteady Mete and tbe bareelkak or Besnby wave 
iMabUlues do m( Intervaae b tha dyaamls. 

In tiia*eelf>|pvrlteting discs, regular and long>llved an* 
tByeknk vertUes have been abeam t o be effldeb at tra^ 

•E 


a dust partlcke b (belr co ps (e g .|Jcb aM ^_g(^li 
eu M slbomdlFii et allODUlgbu et aJJaOldi In 


Sec* 

T? ae (be pertlde (rappbg capablbUs of 

these unsSe^y and Irregular aatkyeksk vortwes. 


9.3 Partlelo ConeontratlOB 

Once t^ gu bjB reached a quasMlsndy gr*/l(o(urbuleBt 
sate, tbe uroJynamlc drag Cute w»d tbe ccsresptsdlng 
hsiVrwTffit taru are btroduced bdo tbe particle and gas 
evduibn, respective^. Ibe u Umb evolved for 


e funber 30 ortdCs, et a*fakS tbie tbe partkka bave bat* 
cally else settlsJ into a quasi saady stale, 'ne parlkls self* 
gravity Is tbeo btroduced and tbe evolved 

far another 6 csfalts Aitnr tbk [lubt, tbe blgb particle denai* 
tks Kbte««d b bound diunpe {sm beloer) r^uR b t^ sdt* 
gravity t«rm in tbe ^bentum 

eqvntbn rnim i|iii'iillj. tbe sep cd tbe tiniulatbb bw 
ceiiMB very sucb that U B Inpractkal to evolve tbe 
partkdes tartbsi Howevts, bspeclton cd tbe den* 

sUy ffMd up to tbB pMnl tbe partUks have 

been suActeotly evdved to draw condual^s about 
tbelr turtbs bebsviouf At this tlioe, particte dynnmKa u 
pc^cUcally no bnger afieSed by gas and one must also bear* 
porate particte sfltsbos. wblU more m 

hlgb particte MtblB ea<b dunqi, to ccsreclly 

falfaw (be subeeguent dy *”**^^ and coauacticn of tbe latter 
tlobanamet sLIImI^ 

surbca density of tbe parttetes at tbe end cd tbe 
run wsb only tbe drag farce <Le, befcselaWo d'<‘"»*rpM(iele 
B pktied In Pig SI As sbows b lOlbb^etaD 
Ild l. tbe oiaD* Md inrersiediste sfaed paitlclei 
wltb UcUoB tlrws f/ a fOOL.01,10(n~' are elBekBtly 
trapped In tbe denSty wave etrtMiitas b t^ gu, wltb 
tbe partkts cd stopping UbM tf s iOfi*' exblbtlng tbe 
bigbest concent rat too Of partkiilar Inlenat far tbe frrimt 
wmIc Is tbe behswur of tbe diM pvtKks b tbe vkln* 
Ry of scale uortves tbnt b aud around (be 

denSty wave strtvtnree. lie dynainlcel t^bavlour d par* 
ticks b tls vtdnlly of «ti— veatkee appenta 

to depsnd eei tbelr stepping tloe, f/. He partkls wltb 
f/ s [0.01,0.l]n~' to Muoulale around tbe cenLral 
of anUcytksb utfUces, but do m( drift farther Into 
cenirea, unih'i'alilr vbda at tbe centres, wblcb 

co n espoi^ to the kcal •*<<»< of PV (^g Q Tbeae parti* 
das b fact (TMe eat tbe over'dense structure around uw« 
UetB, wbkb. es a result, appear to be surrounfad by a spiral, 
or ibg, of dust particks Tlsoe rbga abo lend to 

of tbe higbeel mesMiatlttis of paittrlas b tbe do* 
Bub, (hough bo( slghlflcani^ larger tisn tbe eOBcentiSkss 
faund b (ha ctats of daoulty wavM Thus, tbe —>i^llt* partW 
ds wlib ff= [Cl 01,0.ipi* ^ {weferenilaUy bap out (he vot« 
tlcal strudursB. Ob t^ etber hand, t^ liui iiia ilieli' ilinl 
partitas wiLb Crtcdbi r/ a liO*' get u^pad liBo 
esests of dttsRy waves and evea bwe b Begat Ive PV regMsB 
wltb celallvely low by etoolute value PV (curly arees b Pig 
O, wbkb produce hlgbar over-daslties In tbe gM <Plg (Q 
As opposed to tbe — bicticn nw ease, '*—*» panicka 
811 tbe cestral parts of over-dase r^bos and do bo( kove 
v<^ large partfades wuh Tj = [10 0,100 0)n~‘ acct» 
bulate only b density wave cjtati, v^Uees do not seen (o 
nSect «>!■*» M la dear Irua tbe surface dtsslty of 

t^se particles. 

Hb dlfferwBlal of partlcke b (ha w* 

deal structuras spending on friction tliM can be ex* 
es fallows As doocrlbcd above, ibsue sinaD*scak w* 
dees are iranstent slnkturei, wbkb undergo cyefas of far* 
aaiH^^‘00 B»d tV erlfig away over a fa* orblu 
a> PWi^EBoft Tie siaall perlkks. wttb frk* 
don tlEOes ksB (ban (be orbllel ve dghtly coupled (o 
tbe gM and tberefare cbedy folfaw tbe density vwIsUcbb 
by vortkes. Dtirbg sboA «i*«— tbsoe (ranslenl w* 
dees are appccsdrnalely b by eelf-grsvity, preosure 
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Fi^oM 9. l^^mUmir doai^ erf tlte adfvtdvk) rfiat parade epeOa* {la ibbu of ia itrit {IV) * 9D g b eti Baee i6a 

ba biitLii tbe ^ aad pMdH «wlaBodaced. (via r, b 0 0U)~>) an Aobb latte tep Idl paad aad 

tba odis panek, g to ooeosa^ EnctUo tUaa, toDo* ia ordo wHB tfee laegcet partUhe («iUi t/ m lOtf] ~*) 

la tba bottcea Mt pand tola) etafaM diad^ ol tba dwt >• Aobs la tbe tonuiu QflM ceiapariica, tte m£iM deocr^ 

«f tbe ^ at thM tm a aMaa ta Fi« Q Md the cvT^oadlaj PV 6ehd ■ Aon la F \4 (3) We o^m%« thM 1 m^ paraciM Bttb 
Tf — [10, l(lDfn~‘ are accunwlSed Bi dand^ nee ueaia, Bsmediatooked paitlda wMi b 1 CO~' are c^rUraj afcot»a J y la 
oMd of duty n%w aad ta tte inr|iilntlj iliif il aadcp^mc <.wtk«, vbibt Ae — pMKlM aiib r/ a [(LOi.0 l|fl~‘ —^ w 
teaferatully ~~'*^* aranad aaScpcIoalc forBM bat do aot SU tbau oiml repae, ae laAcatail ^ black vdda, and Baca oot 
i » jj iB [prT 












8 GiA6ou tt of. 



4. UtSmiin Mriaa 4nat7 o( Um j **-'»* nikiii Um 
^ of ume Tb juAip b> tb* dao^ W C a 

107««» uuiMtvsdftte Om mBodcOMa of tbo dra^ forab^naM 
tM (M dual pwttda, vbilft tba scccod /unip ladicMa Iho 
imroductm of tb* ptrucJa Rlf>gtonQ’ 


4bd nXMtta, ^»-«ing ludar'^eAM ntf«, laUb (he 

-iinimum of pv, by flver-deBM ib^ Ae a 

suit, a«n«r pcrtklcB, «4jtiaiKg to ve^tMal 

»»>«■>*« get taddy eo o ca u MW^ Is (hew over^eus rlngt, 
laaotag tbe vonea c«fiir« ^wt^ally eaply. Ob (he o(ber 
Jmiag cehUal (he vcr(K«e svolW, Ihetr PV r» 
dncs by abeohite value aod the Cttrspoadlag over-daitfly 
becBi&a Ar«eg«t abd aa laager hM aa uB^f-deBae core. 
IIm parikkw auh frltlicA Uiaa coBniPaiahfe to the orbital 
>■>■»» 1«nd to acflUBnlate U these owtf-desiilea repaaLlag 
their shape, though with a iKtIe dele?. Larger pertkka are 
Mt afed^ by tbeee traaaleai ■Mt ^ e a . becetwe thefc IrklkA 
Ubb are loa^ than (he orbital (kae - the ehaiactarlstir 
UiiK<srale of vwtiess. 

Figure S shows tla evehuks of the -«a«"riinn. wriue 
density (relative to t^ aeaa) of b the dta 

When the drag (orce batsreefi theparUcIsB andtbe gas 
Is introdiKed after 10 artdU, the surface denbty 

rapidly )imps, lb<Mebg by a of ‘N' 10, boa a tew 
to lOo titaes (be tasea over tlKCOtuas of a faw 

(■him Ilea, the partKks reeoab b a qU0l-ei«ady eiau, 
where they are trapped b vert^es abd/or deaeUy waves, 
bat these ssruidiires gel cobbuouely sbea^ away 

and rsBppenr again over irsr orUtal tbws, so the pertkks 
lesrbg ■ gaseous psrlurbal^B ace swept up by seat esu 
and so on. Thb eyefe coetbiM >'•*>1 the self>|M'Uy cJ the 
part kies Is bt reduced efts dO orUte A( (his point, the par> 
tick concefilraLkbS With large enough tbm> have 

fionoed befcee this (Uk start to undergo rapid ipvrUatlOBal 
collapse, fortabg nany gravKatlonaDy bound doisps of par> 
tides, whkh uulergo turths mergets The laig^ of rk— 
ohjseu exhlUt bcal sutkee densKlb thousands d ttanes 
the BKuo surface desUy These objects are euAdaatly stn> 
bk to ramalo bound after the paceeit denslly wave and/or 
voetlcal perttirbatka they fbmed In had bees ebaared awe?. 
The boundcdnttB of the ■ establubed by eslctilat> 

Ing MtaL hustle plus gravttdlooal energies ^ ||>. 
dividual dunp and, If the latter Is aegatlve, the chuap k 


cksabed ^ pavitatuuslly bound The map of the surfaa 
dei^ty d the at tlM end d the w— i1«iIa»i |g ghown 
B Fig g where we ckaily see bound duni^ of paiildei 
femal due to own sel^grewUy white dots. These 
dunpe are nrostly up ol Intrnnrrllilr tittri panidtn 
with r/ s |0 L, i Ci]f)~‘, w Si)- 4 i tend to conceAlrale Bust efr 
factlvely (see also Fig. |S) laaBes d these gixwipa 
ty^iiaJly d the cedar d 10~^ — 10~’hfs«-es fOlhboae at all 
BSlth AlthBUnM, the gM density and PVdelds an shcFwn 
b Figs (Send Q respectively BeBcarhab^, najcelly d these 
dunpe are (oriMd Mthln the vestkes. with tbe 

above case b^ue seiuhlng ou psrUck selfrgrsvUy (Pig. g), 
t^ pankhs with f/ B 0 If) expedence eJuiaplng in over* 
rings around vertkw, whsreac larges partieks with 
f/ s lJ)d~‘ In nuie uortkes with negative and 

h|p aLsobta value PV (cur^ acees b Pig ( 7 ) As 
ti^<iiir,«j >4 ihowe, these ewdvad voetlMS tand to |fve rise to 
highs «»*-. «ta«iaitUa B ^ tlub t^ ^Bslty woves do. 
B elites frees Fig g (especially the coe In the eaaual* 
lower regloa). Thus, irregular, vsUces appear to be 

Dcee ^benl at trapping parUcks than relatively large 
daoBlly waves. TUa suggests vcellcslparturhalUM can 
be u u^wianl b daoslly wwes b accuasilBing b 
selt-giavllatbg 


4 SVTMMABY AND DISCUSSION 

In this p^HS, we InveBlgstad (hs pBtkk trapping proper* 
Iks cd vorikss and the poteatlal rok they can play b tbe 
(cRuBloo procBS al (he early stage d protepkn* 
e(By dUc Bfe, whn UBl d Ks btarka k »*»■'“ nd csol 
and selPgravliy k a ageeit govecnlag tbe dynunks 

and evolutue d the dlsn Boawvei. bm accretion due to 
t^ grtP/lSatkDal lutabllUy can raise t^ lauperature obd 
activ ate the -'»B»'^*orotatbBal iMtahflitv b tbe »»*>— disc 
(e g, ghjet allSittL Bla . iRicefe kf^i^v<iA but this 
t&c( k bejBBd tbe scope d the present ana^sls, whkh frw 
eitfs ou the eSects d the gratitauoeal iBtablluy cnly We 
used a bmip as Copied b our ptevkus 

WMk (QlbboM but tockasntalkrslBala* 

llou whkh allowed us to capture amai* rrrfilr voetknl 

simeturs against s bBkdrop d laiger scak density waves. 
The laigc seek pkture d the arly stage* d » 

and planet iiw-itinS i«i rhi above papea rematas 

ubchaagod dust p»*eifi— wlihb a quest’Sta^y 

geseons spiral structure produced by a cotohlnad effect d 
diac salf-gravKy and ccolbg The week perivined here gives 
us deeper tnalght into hear stBolhs scak processes aSed 
Um dynacala and evolutbn d paiiklts b the ^u^staady 
gravKoturbuleuce. In this stale, t^ staaller turbulent 
vortka] boUibs (eddls) Issue their traces b the evolution 
d the gas density by ^odiklng ovar-^uss regbos alcng 
with spiral density w a ve s Partkks ore accuniulaled both b 
crests of density wavee and In the antlcyckuk itregulorly* 
shaped uBtBciy vordcB. Vdtkes with sniaD, by aheolute 
value, PV are daiBtarlned by outS'densUlB '*»»•» 

these with larger, by atonlnte vahis, PV and rniininiii'iiillj 
t^ f—eii*ia concent rat tons In tbe feeuMS (urn out (o be ac* 
CDrdmgly higher, even larger lu che daslty waves tor 
r-.aH.. 4ag ifltartMdlato sBod partk^ (Fig g) When tlw 
self.9avtty of the eolld component d t^ dke Is introducsd. 
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FifOM ft. l^^wiihsM suriiM datscr e 4 dtc iadnldiul partieB ipedw (Bi nin* oi A stAA b X^c) m ■ell m (Its WteJ ttvfBM 
(bctMiB o^bC P 4 m 1 ] nt Ute diisi M Um Bid of (M amidMi&fL AssafniES d (Im psads accsdinf to ttoppoig Udm is sobb 
M B F14 S A> ttis RIM, Cbs sv^a dail 9 of ^ s^ (M PV w« tboss, rap^csielr, B Rp. g) and |T] OaJf pw^B wnB 
r/ a ([) 1 , LOjfl upsilBics fvd psvHaOoed pctaiarlty lands iM isspilsr odcTdoBC waUies (ccBpoM triiB Ftp |§sad 

Q. b cnelt vcctsos, psnsciss of (Bs dso to s^sboIms bmi efiRtfsIy Bid tasdi MbttteJr cmimostkIods (sm also Fl$ g), 
'***-^f 10 ihs c4 Bksa* finliiinfiirr pariUs agpaptss, o> (vUis doBJ di—a ^nrviie hn^ aAB 

tM gas uvuidsoBty (b^ b has dlapBBtl 
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6. LogtfiiUiK Mriaa imMCf d iM (<• 4t dte ad ^ 

tbe mil iiUlTO 


pMtkte cpv«r-4eB«IUs fam «1tltlii tbne krr^nl&r 
»««tks Bd «»*« CRMd tapidly uadergo gravltMloul apl> 
l*|Me, bxnlBg *—-iitiip. wttii bnadreda o£ 

Utau laigtr >*»»" ll»e kypl^ gu wubin ^ Aool. 

R«BMtfcnbly, ^ clumps Ccrmed b tbs vert^ss Is 

Isrgn tbsu ihst fcsmsd Is tbs ««v« srasu (Fig 
9 U Is botrswU^ '*»«« tbs dsgss oi psftlds c«s»osMn> 
tioo in tbcac slruciiirts sad bsoes tbsU dumping pi«psrly 
dspsnd on tbal/ frldUm tlam Ws found »>i«i awi-ii aad 
lBl«nimdlats-di 0 d dust carrsspobdlag ta trlclkm 

HfiMi* ccBBpsnbk t«, ot tarn tbnn, tbs local orldtal psrud 
cd tbs 1 */ < 1 (m~‘, bsccims tapsdall) Ugbtly trapped. 
Tlds suggteu tbal MJtlUal perturbstkas, wbkb ennstilly 
are an i wrnitnl ingredteiit, tngrlhrr wlLb denatty avss, la 
t^ (lyoAiBks cd ssK>^avltMlng dlss can also pb^ aa lmp«t> 
i-Ki la tbs ttappbg oi dtM partkks sd plaaslBlBial 
ih^nminti at stages of disc Ills. 

We aould Us to sustn emes again tbal tbs principal ma> 
«r»i tar kigti ds iBi leu -iitiilptig pfrtlplft occairttng la 
tbls pnblsm H tbslrself^avliy Bare, we base not ladu^d 
partkla a ad ona p«i|^i aak boa r tbls slsct esuU 

amdl^ tbs teeulta. (Johaasen w all iOOIlI l bvsatigatsd tbs 
cUert af colbakas B tbs dynamks la gmajue clisa. 

TIk} found tbal oUalaae, wbkii bscssns mms ImporiMit 
al klgb partKls ccaceatratWas, Mtually isad to pRunoia 
partkla cenesnttatlon via dampbg tbalr rta ■sladibK tbs 
^aslUm oi parttfbs bclitdbg eaUuiane liiusd 
out to be lasre than a Cmtas of tbiee blgbar e^parsd lo 
tboss sdtbout TIk i«iUBe<i— dao essm to bavs ao 

^param sSset oa lbs dumping/mllapes ^oesB of 

tbow pmtkk (onesntratfoas dus to tbalr owe aelf-giavtty. 
Tbs masRi o( tbs moel aiasaivs partkie dumps (plaaeUm* 
Mb) kxinsd are talatlvsly basasltns to tbs pseseaes m ak 
iBnoscdodlWcBU Baasd cb tbtee tosulta, ass antklpats that 
Indudcp cd *tJiiBi^« la our ^blem would aol bmA altar 
a bask dyaamleal bsbawltf sd pmtidss, altbougb 
cd Ibis eludy Mtb trsatmant of >-diiain»<a bslmM_Bum^ 
^riipsr>partkka (&> fueb as la lJobanaaa at all 

gOl^ Is oes^d to support tbb conthuBoa. 

la tbU stt^fi ws bavs raasldersd tbs ilmpBdad caas 


Figura 7. n Md at tMww tiiM m ia Fig |5| vfocbi 
ibat la Ftg (Z) Tbs cdwiriiaap la emtnewd agda to Uta naga 
-30 < f < a Tbs ir^M oi tbs PV dnyiiadua we cfoarlr 
aa in tbs —dsBal9 mapa at gsa aad rerTIrlsf ■ tbie tBna 
tepsewd IS Fl^ Sls^m rmpocasdy 


of a rast«tbln (2D) disc wltb a tnoUag Is* la 

order to gain ftrst lasiglit bte tbs sllscta cd vertkss 
on dust partkies la Ibe pits mice cJ selbgtavlly. Obvb 
eusly, far a fulkd undSBlandlag, tbras-dltmaslcnal <SD) 
fn iiiiiiiiil ei vertlcsB, vttb smbsdded partkIsB, la sell* 
grantotiag disee la un i wtij , wbl^, m br as «« are 
awaru, bas not bssa )«1. Tbs dynamks d SD vot> 
Ueas la non<eslfrgr*rllatlng di sa bas bssa iBVSs*‘p’'»*"i 
In a munbtd cd l- 

ibsast all 5553. Lout fc PanaloteoU iMtS Mahsut'eraS 
M12dlP»^>'»T lst alJa)ia hubs 3D ms, a CDolrBA to 


1^ 2D obs, tbsy are subjon lo tbs slBptkal Insublllly, 
wbkb svatnaDy deauoys vwtkss. Besrsvst, stratldcatba 
and barodlale/lisMby vavs drlvlag cm couaterM tbs sl> 
HpUtal iMta^ly, anting vortkss fresn dssuttella A sdt- 
gmvltntiag 3D disc Is clsaily fliatlbsd In tbs vertical dime- 
Uoa. Oa tbs 8elf.gravlty ^«s not fovour lon^ 

Bvsd re^ilar vcrllcai But la tbs SD cass, tic effect of sdt- 
gravi ty ■ irdi>~ibat rsdaesd '^^pnrsd to In tbs 2D cos 

(s g IGoldrsfa-b ^ T.-„4t.wBelJlui55lM»*««»^*'ttH gHa 

i20lffl Tbus, K Rtnnbs tobssssalafutursruuarkaleludleB 
bow Ibe eocoblnatmo of tbs sbevs foclom ebnpes tbs dynnea* 
ks of SD vrUcss b edf-gravitaUng duta Wbea aanlytbg 
bshavltjur cd parudes b audi 3D wtkM b aelf>gravltailng 
dlscB, oae must also talw bio account tbs bask eSreU b 
1^ particle dynamks at work la SD - ti iliiia iililkm 
dus to vertical giavit}' aad ellr np and dumping id tbs dad 
layer dim to tbs Kdvla-Hebnbdls and slnonlB^BS^U- 
tbs, fsepscMvsIy (s^.'^udiai-SbillSBu**-****' ■» 


Jobaoea le BOOT. ai jilJ IjpOa 


Tbs wwk |iii'eriiliii 1 ben on tbs pldun of 

lbs early slagta of plaa *-^*--' forroatfoa boib up by 
labbtas el air Bildl . wblcb dtlD pteasald aattn» 
Uvs msUmd for lbs rap^ crsatlcn of a large leservMr af 
r*—afoag altb arvaral iiiiieslirn objects «f man 
‘N' 10~*Afs«Fft Tbsae ob)ec1s are Hkdy to 
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to gfc— lia tbt tt »^iihr>«nJ eot« AccntlCD praccat (e g., 
Rmber •wk. **Uekaitais (b« 


iPolhtk M J 1 T 5 

nsat of Urg», gtniUnll^ofirj bound ^ 

tides Mik 'eiAk' pwtkiee, InlmcUag aub ibe disc sad 
arlih each other, ooiiM expand the scope ^ gniiiU»b>«ifl flu-fc 
M thcas proseated ^ovn This would sllow £■ U»e sub* 
■a»'» ovoluUu of these plnnetesual-elsed objects nod 
could {soved hirlhts kalgbt bo* qui^ly the cote a» 
crotiOB procBse can oork In aelf-gi6vlt«lBig; 
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APPENTklX A» RESOLUTION TEST 

Tb lBTCeU£M« Ibe «8ecu tt naoliiUiei eu (be dyaeaitfs of 
teruces In the piumec cf self-^avU}', we teu g<rnniib>«ia 
So€ ft denftlo wHb ball Mfte, L,/i = ^/3 s 10/« (la 
ct H. Biltlftl^ B1) ftl three dlftareat 
266x2S6,Sl2vSl2,1034 1034, bcludlog the Mine umber 
c4 perUcia and Mtb tbe inlUnl 

M Ui Ibe pftper. Ic equlvekut t« raanlnf tbe Odndftl 
Bodd wKb La s Ly s 30/« cAiklerwJ In LbU work eC 
doubkd laftlntkBft 612 s 312,1024 x 1024 end 204S x 2048, 
reepecttuely, tbe Mnaponding 1^ ceD fieee. 4s s 
La/N, B 0J)134,0.0062. (L0031 weuU be tbe nm In Imb 
rng~i Se Car a umerlcel cmveigesee 1» ccaeemed, gnd 
ceD sine It a eential leb^h>«eftk, utBCub^ warloue 
pbjftlenl len^lis df system Kgaixai K sUaws us te 
cut wbetbsi a aBnulnUoa U resolved '^erefore. diK to tbe 
ptd (ell stMS, a tesolutkn stndy we do fw tbe balLefete 
m>del also atabUsbs wbetbar vortlcs In tbe dduccal 
are suflciently l asal^eA 

la tbe qUBSi-eteady state, tbe PV Selds tnm aO tbese 
three runs at tbe ante tone are sbowa In Pig. AL Tbsa 
fleldi, at dlSecent rseoktsoBB, dl^pls^ sasenLlaDy 

Identical structures, '*««» tbey are resolved 1b 

PBi ^ m ue quanUtallve grounds, m in 

IJo b a sn n Ir Gammlfl 620061 . we consIdB’ tbe autoctstelatloo 
AiDCtbn tfA tbe ductuatbg part o( PV, df = I — {/}, given 
hy 

C(A*) = (#/(*)«/(»♦ A*)), 


wfanre tbe "xc** kackeu denote average ever tbe whole dtv 
-««■>« Figure A2 abawa tbe antoeorrelatloe tuncttons cor* 
respondkag to tbe PV HMt In Fl^ AL Tbey do appear 
■Sfiiiu- and tbeir cors bave nearly wi alUptlml shape with 
aflnu and nake principal ■*»« FoUowtag tbe of 

lOnan at J1120090. we calculated related ees’ i e l atlon 
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Figure Al. FV bid to tbs tell tte tesnaln at difesB 
Boos in tbe qBaw>4tmi^ state at tbe sama dma (cdoor aalaa) 


tbsve aus by 

“ mC 

wbre i ■ tbe frimi Ax s 0 tbe principal axM 

given by tbe ‘mii ewclcn Amv* and Tlse OAelatlon 
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Isaftlu CAb b« i««d at »laoiBiiK fix tiM (ypk&l aiM of (be 
PV sLrucluna 

Tbe soelleet correlftUoa kngUia almg tiM uli 
are W » 0.088. (LOn.OOn, fix (be lo*W. UU^oMdlftte 
and nsahiUaB rune, respective, 'ne ratloe ei (be 

-^mtnntn knjtb to »*i» grid cell ■( (bese 

rsahitaans are, rspeeikvelf, Xm^/Ss = 7 1,1±7,24 6. Tbk 
nLla boeases eHtb resdhitdCD wblcb serves m an 
fiM cttvergeace and, as a mallet <d Wt, Implies that (be 
voetkes (eddks) ua be maldered to be suffidjeady r^olved 
lb (be tnodd oblA bee Sx = 0.0062 In our {xevkpus 

fteaeea lGihbe^etall <511 , BH^ . tiM domain abe Bl. s 
L, s 8D/e abd reaDbelMb 1034 ^ 1CB4 fix (be earns o(be 
panmeleta, giving it = 0.0240 and bean Xmim/it < 
Tbue, altbougb IbB resohiuon is q«d(e enitable to resolve 
large scak ^nslty emvea b (bme eludhs, K la sUU coarse lo 
detect and exwnlne vortical eiructurea at wkL laige 
eUs b ccDUast lo (be al four Ui^ bcJW gua 

used b (be picerni pape. 


'/«. ^ »6 



41S -91 -OH 0 OK 01 0.1* 


M 



N. 



N. 


^gure A3. Nenulbed autoaereBoo fnacOnn, • vbcre 

— <(0), to (he PV la Pig. Ai M dtSxaB 

reaobiBMia feclots enhae) 
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